Quantifying the intrinsic amount of fabrication disorder in photonic-crystal 
waveguides from optical far-field intensity measurements 
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Residual disorder due to fabrication imperfections has important impact in nanophotonics where it 
may degrade device performance by increasing radiation loss or spontaneously trap light by Anderson 
localization. We propose and demonstrate experimentally a method of quantifying the intrinsic 
amount of disorder in state-of-the-art photonic-crystal waveguides from far-field measurements of 
the Anderson- localized modes. This is achieved by comparing the spectral range that Anderson 
localization is observed to numerical simulations and the method offers sensitivity down to c± 1 nm. 



Optical nanostructures are able to control the emission 
and propagation of light at the nanoscale.^ Such structures 
are generally made by careful engineering where fabrication 
imperfections give rise to ubiquitous disorder, which may 
degrade the performance of nanophotonics devices inducing 
losseP^or Anderson localization of light Given the matu- 
rity of current fabrication processes, quantifying the intrinsic 
amount of disorder in photonic nanostructures is challeng- 
ing since any imaging technique would require subnanome- 
ter spatial resolution to accurately resolve the structural fea- 
tures. Important progress has been made in order to mea- 
sure such unavoidable disorder by, e.g, statistical analysis of 
scanning electron micrographs. 8 However, reconstruction of 
actual dielectric profiles from scanning electron microscopy is 
not straightforward. 9 Here we demonstrate a way of quantify- 
ing the intrinsic disorder of photonic-crystal waveguides from 
optical far-field intensity measurements. Our approach offers 
an accurate statistical measurement of the imperfections with 
nanometer resolution. 

As initially explored in the context of electronic transport, 
the general eigenstates of a perturbed Hamiltonian describing 
a crystal affected by structural disorder was obtained by I. 
M. Lifshitz. 10 The energy spectrum contains states that are 
different from propagating Bloch solutions and, on average, 
decay exponentially with a characteristic length scale called 
the localization length. These impurity levels populate the 
spectrum near the band edge of the unperturbed crystal with 
narrow resonances and smear out the density of states (DOS) 
by forming the so-called Lifshitz tail. This situation applies 
not only to the electronic wave case but to any kind of wave 
propagating through a perturbed periodic lattice. The optical 
Lifshitz tail has been probed by scanning near-field optical 
microscopy^ and by single quantum-emitter spectroscopy™ 
The width of the Lifshitz tail is determined by the impurity 
density 10 or, in nanophotonic structures, by the amount of 
structural disorder. This dependence can be used to quantify 
the residual disorder of photonic nanostructures. 

We use two-dimensional finite- difference time-domain 
(FDTD) simulations 13 to calculate the Lifshitz tail of 
photonic-crystal waveguides using a freely available software 
package™ The simulated structures replicate the fabricated 
ones and consist of a hexagonal lattice of air holes forming 
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Figure 1: (Color online) (a) Two-dimensional finite- 
difference time-domain calculation of the electromagnetic field 
component perpendicular to a photonic-crystal waveguide 
(color scale). The parameters used for the calculation are 
given in the text. Disorder is introduced in the structure by 
randomly displacing the holes position in the three rows on 
each side of the waveguide following a normal distribution 
with standard deviation of a = 0.04a (see inset), (b) Disper- 
sion relation of an even-parity photonic crystal guided mode 
(dashed-red curve). The shaded region around the mode cutoff 
outlines the spectral region where Anderson-localized modes 
appear determining the Lifshitz tail, (c) Calculated normal- 
ized electromagnetic field intensity ensemble-averaged over 20 
different positions along a single waveguide and a total of 10 
disordered photonic-crystal waveguides with a = 0.04a. The 
residual speckles in the spectrum are due to the finite statis- 
tics. The spectral width of the region where localized modes 
appear can be determined precisely and used to extract the 
degree of disorder in experimental samples. The inset shows 
the Lifshitz tails for different amounts of disorder. 
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a two-dimensional photonic crystal with a lattice constant 
a — 260 nm and a hole radius 0.29a. A waveguide is intro- 
duced in this structure by leaving out a row of holes. The 
length of the simulation domain is 120a, with eight rows of 
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Figure 2: (Color online) (a) Scanning-electron micrograph 
of a photonic-crystal waveguide (top view) without engi- 
neered disorder. The inset illustrates that unavoidable irreg- 
ular shapes and sizes of the holes are present in the samples, 
(b) Normalized high-power photoluminescence spectra col- 
lected while scanning a microscope objective along a photonic- 
crystal waveguide with only intrinsic disorder (cr e = 0). The 
bright peaks that appear randomly along the waveguide and 
close to the waveguide mode cutoff correspond to Anderson- 
localized modes. In order to display the peaks of varying 
heights, the plot is saturated at 10 % of the maximum peak 
intensity. The spectral region where the modes appear de- 
fines the experimental Lifshitz tail, (c) Histogram of the spec- 
tral position of the Anderson-localized modes for different 
amounts of engineered disorder, a e . The width of the Lifshitz 
tail increases when increasing the amount of disorder. 



holes on each side of the waveguide. The position of the holes 
in the three rows on both sides of the waveguide are randomly 
displaced a certain amount, Ar, normally distributed with 
a standard deviation a = V< Ar 2 > — < Ar > 2 , where the 
brackets indicate the expectation value. By using an effective 
refractive index of n = 2.76, the 150 \im thick photonic-crystal 
membrane can effectively be simulated in 2D, 15 -' which signifi- 
cantly improves computation times. Perfectly matching layers 
(PML) cover 10 lattice constants at both waveguide termina- 
tions in order to model an open system.^ A dramatic effect 
of disorder in a photonic-crystal waveguide is one-dimensional 
Anderson localization, which confines light randomly along 
the waveguide forming a narrow spectral band of strongl y lo- 
calized resonances around the waveguide cutoff frequencyP^ 
We use a total of nine electric dipole sources randomly po- 
sitioned along the waveguide to excite the localized modes 
in a broad bandwidth of Acj = 0.02(a/A), centered at the 
waveguide mode cutoff frequency. From the resulting elec- 
tric field, E(r,t), we can calculate the intensity spectrum, 
= | J r J r T[E(r, £)] | 2 that consists of a set of Lorentzian 



peaks due to the Anderson-localized quasimodes. To circum- 
vent the limited spectral resolution of the Fourier transform, 
we use a method denoted harmonic inversion, which allows 
resolving spectrally-close resonances and to extract their Q- 
factors with high accuracy.^ Further details of the simula- 
tions can be found in Ref. [18]. Figure TJa) plots the in-plane 
component of the electric field perpendicular to the photonic- 
crystal waveguide, E y (r), revealing strongly- confined modes 
randomly appearing along the waveguide. 

Figure [TJb) shows the dispersion relation, u(k), of the 
unperturbed structure,^ where k is the length of the wave 
vector along the waveguide. The fundamental TE-guided 
mode (dashed-red line) gives rise to a strongly increasing 
DOS (DOS = (l/ir)dk/dcj)j enabling slow light with very 
low group velocity (v g oc 1/D OSp^ and efficient broad- 
band single-photon sources fi 2 * 21 * Residual disorder smears the 
DOS 22 , sets a limitation to slow light, and may induce lo- 
calized modes in the slow-light regime if the sample length 
and the characteristic length related to out-of-plane losses ex- 
ceed the localization lengthP^' Figure [ljc) shows the calcu- 
lated ensemble- averaged intensity, displaying the spec- 
tral range of the Lifshitz tail. The observed sharp resonances 
illustrate that the finite statistics implies that the ensemble 
average is not fully converged, but nonetheless the width of 
the Lifshitz tail can be precisely extracted. We define the 
width of the Lifshitz tail, Acj, as the spectral range where 
sharp resonances appear, i.e., peaks that are, at least, 10 times 
brighter than at neighboring frequencies. As shown in the in- 
set of Fig. [ljc), the width of the Lifshitz tail broadens when 
increasing a and we can use this dependence to extract the in- 
trinsic amount of disorder in the photonic-crystal waveguides. 

Lifshitz tails are measured in 100 urn long GaAs photonic- 
crystal waveguide membranes containing a single layer of self- 
assembled InAs quantum dots (QDs) in the center and oth- 
erwise identical parameters as the simulated structures (see 
Fig.[2ja), details of the samples can be found in Ref. [6])- Due 
to the tolerance of the fabrication process, imperfections in 
the shape, size, and position of the holes are unavoidable, see 
inset of Fig. |2ja). A controlled degree of disorder in the hole 
position is introduced in the three rows of holes on both sides 
of the waveguide in the same way as discussed for the simu- 
lations. The embedded QDs are used to excite the Anderson- 
localized modes by using a confocal micro-photoluminescence 
setup for excitation and collection [211124] . Photoluminescence 
emission spectra are collected in a wide wavelength range 
of 900 nm < A < 1000 nm while scanning the microscope 
objective along the photonic-crystal waveguides, as shown 
in Fig. |5Jb). Clear Anderson- localized modes are observed 
within a narrow spectral range around the waveguide cutoff 
(A = 978 nm), which provides the width of the Lifshitz tail. 
Figure [2|c) displays the histogram of the Anderson- localized 
modes spectral positions for different amounts of engineered 
disorder, showing how the Lifshitz tail broadens spectrally 
when increasing the amount of disorder as expected from the 
simulations. 

For small amounts of disorder, the calculated width of the 
Lifshitz tail depends linearly on the amount of disorder, i.e., 
AcJcai oc a, as displayed in Fig. [3^. Note that in the limit 
of very small amount of disorder the applied boundary con- 
ditions do not suppress stray reflections at the end of the 
photonic-crystal waveguide sufficiently well, which appears as 
a finite width of the Lifshitz tail. As predicted by the simula- 
tions, the width of the Lifshitz tail recorded experimentally, 
Acj e xp, is found to vary linearly with the engineered disorder, 
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Figure 3: (Color online) (a) Calculated width of the Lifshitz 
tail for photonic-crystal waveguides with a different standard 
deviation in the hole positions and a — 260 nm. The fit to a 
linear dependence on disorder (red line) is very satisfactory 
for small amounts of disorder (< a — 0.1a). (b) Experimental 
width of the Lifshitz tail measured for different amounts of 
engineered disorder, (c) Recorded intrinsic degree of disorder 
for samples with different engineered disorder. 



see Fig. [3]o . It turns out that the Anderson-localized modes 
appearing close to the waveguide terminations are spectrally 
shifted. In order to avoid such finite-size effects, the width 
of the Lifshitz tail is calculated by dividing the waveguide 
into sections of 20 \xm. length and integrate the spectrum over 
each of these sections. The mean value and standard devia- 
tion of the Lifshitz tail width are shown in Fig. [3b for different 



amounts of disorder. These data can be used to extract the in- 
trinsic disorder associated with the fabrication process. Thus, 
the total degree of disorder in the samples a consists of both 
intrinsic (<Ji) and engineered disorder (a e ) that are assumed 
to be independent and therefore related through 



2 . 2 
(Ji +CT e . 



(1) 



Making use of the predicted linear relation between Acj and 
a enables us therefore to extract <ii from the data in Fig. [3J). 
It is displayed in Fig. [§}; for the various samples with different 
amounts of engineered disorder, and as anticipated no corre- 
lation between the two is observed. While the method does 
not identify the nature of the intrinsic disorder, we are able 
to extract quantitative numbers for the statistical properties 
of the disorder. Thus the average degree of intrinsic disorder 
found in the samples is equivalent to hole displacements of 
< (Ji >= (1.2±1.0)nm. These results appear to be comparable 
to the numbers inferred from optical near-field measurements 
of the Lifshitz taiP^ and from statistical analysis of scanning 
electron micrographs 8 . We emphasize the remarkable sensi- 
tivity of the applied method since it allows to probe disorder 
at the sub-nanometer length scale. Furthermore, using the 
QDs as internal emitters for exciting the Anderson-localized 
modes appears a very attractive way of quantifying disorder 
compared to tedious statistical analysis or complicated phase- 
sensitive optical near- field measurements. 

In conclusion, we have measured the spectral width of 
the Lifshitz tail of photonic-crystal waveguides with different 
amounts of engineered disorder in the regime of ID Anderson 
localization. By comparing the experimental data to numeri- 
cal simulations, we are able to quantify the intrinsic disorder 
of the samples with nanometer resolution. The method pro- 
posed here could potentially be implemented in other disor- 
dered nanophotonic structures as well in order to quantify and 
potentially optimize the degree of disorder within the settin g 
of cavity quantum-electrodynamics 2 ^ or random lasin g 26 * 27 *. 
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